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Abstract 

By means of the density matrix renormalization group (DMRG) method, the static spin structure 
factors and the magnetization plateaus of the trimerized Heisenberg ferromagnet-ferromagnet- 
antiferromagnet and antiferromagnet-antiferromagnet-ferromagnet spin chains in the presence of a 
magnetic field are elaborately studied. It is found that in the plateau states, the static structure 
factor with three peaks does not vary with the external magnetic field as well as the exchange 
couplings; the spin correlation function behaves as a perfect sequence and has a simple relation 
with the magnetization per site. An approximate wave function for the plateau states is proposed, 
and a picture based on the valence-bond-solid states is presented in order to understand the origin 
and the total number of the magnetization plateaus, which are shown to be in agreement with the 
DMRG results. 

PACS numbers: 75.10.Jm, 75.40. Cx 
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In recent years, owing to their unusual magnetic properties, the low-dimensional quantum 
spin systems have been extensively investigated both theoretically and experimentally (see 
e.g. Ref. jj). Among others, the phenomenon of topological quantization of magnetization 



has received particular attention (e.g. 



jmenon 01 topological quantizat 



13 



14 



15 



16|,H). 



Recently, the magnetic and thermodynamic properties of the ferromagnet-ferromagnet- 
antiferromagnet (F-F-AF) and antiferromagnet-antiferromagnet-ferromagnet (AF-AF-F) 
trimerized Heisenberg spin chains have been systematically explored fl^. fig. For both sys- 
tems, it have been observed that the topological quantization condition n(S — m) = integer, 
with n the period of the ground state determined by the explicit spatial structure of Hamil- 
tonian, S the spin, and m the magnetization per site, is only a necessary, but not a sufficient, 
condition for the appearance of magnetization plateaus, and the total number of the plateaus 
is 25* + 1. The width of plateaus is found to decrease with increasing the ferromagnetic in- 
teraction, and the plateau-non-plateau transition is of Kosterlitz-Thouless type, while the 
position of the plateaus is independent of the magnitude of the exchange couplings jjjj]. The 
magnetization plateau is smeared out when temperature is increased, and it appears the sys- 
tems fall into different thermodynamic states under different magnetic fields[15(]. In order 
for understanding the spin-spin correlations in such intriguing systems, and demonstrating 
why the total number of the plateaus is 2S + 1, we shall, in this paper, investigate the 
static spin structure factor that would be directly measured experimentally for F-F-AF and 
AF-AF-F trimerized Heisenberg chains with S = 1/2 and 1, and offer a picture based on the 
valence-bond-solid (VBS) states for further understanding the origin and the total number 
of the plateaus. 

The Hamiltonian of the trimerized Heisenberg spin chain in a magnetic field is given by 



H — ^(JS3j_2 ■ S3j_i + JS3j_i ■ S3./ 
3 

+ J'S 3 j • S 3 j + i) — Sj, 



(1) 



where J, J' > denote the AF coupling and < the F coupling, and h is the external 
magnetic field. We take g\is — 1 for simplicity. Two cases with 5 = 1/2 and 1 will be 
considered: (a) J = — Jp < 0, J' = Jaf > 0, the system (F-F-AF) is an antiferromagnet; (b) 
J = Jaf > 0, J' = —Jf < 0, the system (AF-AF-F) is a ferrimagnet. The static structure 
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factor S(q) is the Fourier transform of the spin correlation function (S^S^,), defined by 



S(q) = ^E e^ n - n '\S z n S z n ,). (2) 



n,n' 



The density matrix renormalization group (DMRG) method Jjj, |20j will be invoked. In the 
following calculations, the length of the chain is taken as iV = 60, and the number of the 
optimal states is kept as M = 60. Open boundary conditions are adopted. The truncation 
error is less than 10~ 3 in all calculations. 

For the spin-1/2 F-F-AF chain with couplings Jaf/Jf = 1, one magnetization plateau 
with m = 1/6 is observed, being consistent with the quantization condition [5[]. Let h c i and 
h c2 be the lower and upper critical fields at which the magnetization plateau appears and 
disappears, respectively, and h s be the saturation field, as indicated in Fig. [T](a). When 
h < h c i, as demonstrated in Fig. [T](b), the static spin structure factor S(q) shows three 
peaks at q = 7r/3, 5tt/3, and 7r. With increasing the magnetic field, the peak at q = tt 
disappears while the peak at q = 7r/3 (57r/3) splits into two small peaks that continuously 
shift towards q = and 2it/3 (q = and Air/ 3), respectively. Meanwhile, new peaks start 
to appear in q = 0, 27r/3 and 47r/3. When h c \ < h < h c2) where the ground states are in the 
plateau state with m = 1/6, as shown in Fig. [Etc), S(q) with three peaks at q — 0, 27r/3 
and 4ir/3 does not vary with the magnetic field, which corresponds to the long-range order 
and a perfect sequence for the spin correlation function (SqSj), which we will discuss later. 
When h C 2 < h < h s , the peaks at q = 2tt/3 and Att/3 are suppressed with the increasing 
field and disappear at the saturation field, while the peak at q = increases with the field, 
which show the ferromagnetic characteristic, as demonstrated in Fig. [H(d). 

In the absence of the external field, we find that the static structure factor S(q) could be 
fitted by 



s (q) = XX a ' e 



i=i 



_ip.2 sin 3g(cosh3/5 — cos3g) sin/g 
cosh(6/3) — cos(6g) 



(e 6/3 — cos 6q — 2 cos 3q sinh 3/3) cos Iq 
cosh(6/3) — cos(6g) 

+ ^til, (3) 

where ati and f3 are parameters. As presented in Fig. [T](b), the characteristic peaks can be 
well fitted by Eq. ([3]), with only a slightly quantitative derivation, showing that the main 
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features of S(q) of this model can be captured by a linear superposition of three modes, 
which is closely related to the low-lying excitations of these trimerized spin systems 

In the plateau states, S(q) has three peaks at q = 0, 2n/3 and 4tt/3, and (SqSj) shows 
a perfect sequence with a period of 3. It is the period of (S^S*) that determines the peak 
positions of S(q). Define £(m) = (l/N)^2 n (S^S^_ m ), which also shows a perfect sequence 
with a period of 3, say, {• • • , (£(3j),£(3j + l),£(3j + 2)), • • •}. Thus, the structure factor 
S(q) in the plateau states can be written as 

S(q) = ls(S + l) + 2f(l) cos(g) + 2f(2) cos(2g) 

OO OO 

+ 2 £ CM cos(3jg) + 2 £ £(3j + 1) cos(3j + l)q 

3=1 3=1 
oo 

+ 2^(3j + 2)cos(3j + 2)g. (4) 

3=1 

Because the series YfjLi cos(3j(7) only diverges at q = 0, 27r/3 and 47r/3, and oscillates at 
other q, the peak positions of S(q) appear only at q — 0, 27r/3 and 47r/3, consistent with our 
numerical calculations. This argumentation may be generalized to other spin chains, say, 
if (SqSj) behaves as a perfect sequence with a period of n, S(q) would only have peaks at 
q — 0, 2Tr/n, 4Tr/n, ■ ■ -. 

Besides the sequence structure of the spin correlation function, it is found that, from 
our DMRG results, in the plateau states the spin correlation functions satisfy the following 
relation: 

|[e(3j) + £(3j + 1) + £(3j + 2)] = m\ (5) 

which does not alter with the couplings Jaf/Jf as well as the magnitude of the spin S. As 
the susceptibility xCO = k^fJf Sij^f^i) — ( -S 1 ? ) ( S 1 ? ) ] bears a finite value at T — > in the 
plateau states because of m(T) being a constant, it leads to 

^Y,[(S;Sfi-(S?>(S3)] = (6) 
h3 

at T — >• 0. Since (SqS*) shows a perfect sequence in the plateau states, we can get Eq. (jHJ) 
directly from Eq. (jBJ). 

The spin-1 F-F-AF Heisenberg chain with the couplings Jaf/Jf = 1 is also studied. 
The two plateaus at m = 1/3 and 2/3 are observed. The static structure factor and spin 
correlation function show the same qualitative characteristics as those of the spin- 1/2 chain. 
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FIG. 1: (Color online) For a spin-1/2 F-F-AF trimerized Heisenberg chain with couplings Jaf/Jf 
I . (a) the magnetization per site m as a function of the external magnetic field h in the ground 
states; and the static structure factor S(q) in the ground states under the external field (b) h < h c ±, 
(c) h c i < h < h C 2, and (d) h C 2 < h < h s . 

The zero-field static structure factor S(q) can also be well fitted by Eq. ([3]), and in both 
plateaus, Eq. dSJ) is still satisfied. 

The AF-AF-F trimerized Heisenberg chains with spin 1/2 and 1 are investigated as well. 
For the spin-1 AF-AF-F Heisenberg chain with couplings Jaf/Jf = 1, two plateaus at 
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m = 1/3 and m = 2/3 are observed. Let h c \ denote the critical field of disappearance 
of m = 1/3 plateau, h c2 and h c ^ be the critical fields of appearance and disappearance of 
m = 2/3 plateau, respectively, and h s be the saturation field, as marked in Fig. [2(a). When 
no magnetic field is applied, S(q) shows three peaks at q = 0, 27r/3 and 47r/3, which differs 
from the case of the F-F-AF chain because the latter is an antiferromagnet. When h < h c \ 
and h c2 < h < h c3 , where the ground states are in m = 1/3 and 2/3 plateau states, S(q) 
with three peaks at q = 0, 27r/3 and An/ 3 does not vary with the external field, both of 
which show a perfect sequence with a period of 3 for (S^Sp, as shown in Fig. [2(b). When 
h C 3 < h < h s , like the case of the F-F-AF chain, the peaks at q = 27r/3 and 47r/3 are 
depressed with increasing the field, and disappear at the saturation field, while the peak at 
q = increases with the field, as shown in Fig. [2(d). 

The spin-1/2 AF-AF-F Heisenberg chain with couplings Jaf/Jf = 1 is also explored. 
The features similar to those of S — 1 are seen, and the spin correlation functions also 
satisfy Eq. (jSJ) in the plateau states. 

The plateau states for the F-F-AF and AF-AF-F trimerized Heisenberg chains with spin 
S can be described by an approximate wave function defined by 

Ipi = a | i3i-2|3i-l |3«) + b\ T 3i— 2 1 3i — 1 T 3i ) 

+ c\ Isi-sTsi-itsi), (i = !,■■■, N/3) (7) 

where |j (ij) denotes spin up (down) on site j. It is found that, if the coefficients a, b and 
c are properly chosen, one could find the so-calculated local magnetic moment (Sp as well 
as the spin correlation functions to be in agreement with the DMRG results. Take the spin- 
1/2 F-F-AF trimerized Heisenberg chain for example. In the plateau state with m = 1/6, 
the DMRG results show that (S?) behaves as {• • • , (Si, S2, S2), • • •} with S\ = 0.367 and 
5*2 = 0.067, giving rise to the magnetization per site m = (Si + 2S 2 )/3 = 1/6, and the spin 
correlation function £(j) behaves such that {• • • , (£,(3j),£(3j+l),£ i (3j+2)), ■ ■ ■} with £(3j) = 
0.048 and£(3j + l) = f (3j + 2) = 0.018, satisfying I[£(3j)+£(3j + 1)+£(3j + 2)] = (|) 2 . If we 
choose the coefficients a = ±^13/30, b = ±^/l3/30 and c = ±y2/15, from the approximate 
wave function we can obtain (4 , i\S^ i _ 2 \4'i) = 11/30, (V'il^-ilV'i) = = 1/15, and 

£(3j) = 43/900, £(3j + 1) = £(3j + 2) = 4/225, both of which are quite in agreement with 
our DMRG results. 

Let us now invoke the VBS picture to show why the total number of the magnetization 
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plateaus in these trimerized systems is 25" + 1. As shown in Fig. [31 we denote a spin 1/2 
by a bullet, a spin singlet of two spin 1/2 by a short line between two bullets, and the 
spin 5 per site by a large open ellipse that could be viewed as a combination of 25 spin 
1/2. For the F-F-AF chains with spin 5 = 1/2, 1, 3/2 and 2, the isolated spins 5 that are 
denoted by large open ellipses, are separated by spin singlet states. Thus, these isolated 
spins could be regarded as a classic spin chain with antiferromagnetic couplings, whose 
magnetization curves should have no plateau except for the saturation plateau of m = 5/3. 
This picture explains why the magnetization plateau of m < 5/3, which is permitted by 
the quantization condition [a], does not appear. For the ferrimagnetic AF-AF-F chains with 
5 = 1/2, 1, 3/2 and 2, the isolated spins are aligned parallel, and the magnetization per 
site in the absence of the magnetic field is 5/3, leading to the plateaus of m < 5/3 never 
exist. With increasing the magnetic field, the spin singlet bonds are gradually broken. When 
a spin singlet bond per three sites is simultaneously broken, a new magnetization plateau 
appears, and so on. When 25 spin singlet bonds per three sites are simultaneously broken, 
there are 25 plateaus occurring. Considering that there is also a saturation plateau, the 
magnetization per site obeys m = (5 + x)/3 (x = 0, 1, • • • , 25). Therefore, the total number 



of magnetization plateaus is 25 + 1, which is well in agreement with the DMRG results [141] . 
For instance, the trimerized Heisenberg chain with 5 = 2, according to our VBS picture, 
exhibits magnetization plateaus at m — 2/3, 1, 4/3, 5/3 and 2, which have been confirmed 
by the DMRG results 14]. 

In summary, we have studied the static structure factor 5(g) and the plateau states of 
the F-F-AF and AF-AF-F trimerized Heisenberg chains. It has been found that in plateau 
states, 5(g) with three peaks at q = 0, 27r/3 and 47r/3 does not vary with the magnetic field 
as well as the exchange couplings Jaf/Jf, and the peak positions are determined by the 
period of the spin correlation function. From the DMRG data, it is observed that the spin 
correlation function satisfies a simple relation with the magnetization per site [Eq.fjS])]. An 
approximate trimerized state is presented to describe the physical properties in the plateau 
states, and the results obtained by using this trial wave function are quite in agreement with 
the DMRG data. To explain the origin and the total number of the plateaus, a picture in 
terms of the VBS states is offered. From this picture, the appearance of the magnetization 
plateaus can be attributed to the breaking of the spin singlet states and the isolated spins. 
The total number and the positions of the plateaus deduced from our VBS picture are 
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exactly the same with the DMRG results. It is expected that the properties of the static 
structure factors reported in this paper could be useful for future experimental studies on 
the trimerized Heisenberg spin chain compounds. 
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FIG. 2: (Color online) For a spin-1 AF-AF-F trimerized Heisenberg chain with fixed couplings 
Jaf/Jf = 1, ( a ) the magnetization per site m as a function of the external magnetic field h in the 
ground states; and the static structure factor S(q) in the ground states with the external field (b) 
h < h c \ and h C 2 < h < h c s, (c) h c \ <h< h C 2, (d) h c ^ < h < h s . 
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FIG. 3: (Color online) Sketch of the VBS picture of the F-F-AF and AF-AF-F trimerized Heisen- 
berg chains with spin (a) S = 1/2, (b) 1, (c) 3/2, and (d) 2, respectively, f denotes spin up, and 
I spin down. The bullet • denotes a spin 1/2, and the short line between bullets denotes the spin 
singlet state. The big open elipse presents spin S per site which could be viewed as a combination 
of 2S spin 1/2. 



11 



